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Abstract 

This paper analyzes the influence of general, small volume, inclusions on the 
trace at the domain's boundary of the solution to elliptic equations of the form 
V • D'^Vu'" = or (—A + q^)u'^ = with prescribed Neumann conditions. The 
theory is well-known when the constitutive parameters in the elliptic equation 
assume the values of different and smooth functions in the background and inside 
the inclusions. We generalize the results to the case of arbitrary, and thus possibly 
rapid, fluctuations of the parameters inside the inclusion and obtain expansions 
of the trace of the solution at the domain's boundary up to an order e^'^, where 
d is dimension and e is the diameter of the inclusion. We construct inclusions 
whose leading influence is of order at most e'^'^^ rather than the expected e*^. We 
also compare the expansions for the diffusion and Helmholtz equation and their 
relationship via the classical Liouville change of variables. 

keywords: Small volume inclusions, asymptotic expansions, diffusion equation, Helmholtz 
equation, inverse problems. 

AMS: 35J15, 35B40, 65R20, 35R30. 

1 Introduction 



Asymptotic expansions for the influence of small volume inclusions for elliptic and other 
equations is now well-established. We refer the reader to e.g. [21 [3l HI [5l [6l [8] and 
their references for a few historic and recent works on the subject. A major advantage 
of such expansions is that they help us understand what details of the constitutive 
parameters in the equation may or may not be reconstructed from available boundary 
measurements. Indeed, in the elliptic equations of interest in this paper, namely the 
diffusion or conductivity equation and the Helmholtz equation, the reconstruction of the 
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constitutive parameters X from knowledge of the full Dirichlet-to-Neumann map A, the 
most general type of information available at the domain's boundary, is an extremely ill- 
conditioned problem. Available stability estimates for both types of equations predict 
that the accuracy in the reconstruction is at best logarithmic in the accuracy of the 
measurements. More precisely, we have [H [TO] 

-s 

for some positive constant C and 6 G (0, 1), where Xi and X2 are two sets of parameters 
and Ai and A2 their corresponding measurements. 

For such severely ill-posed problems, only a limited number of degrees of freedom may 
be reconstructed from even quite accurate measurements. A natural way of limiting the 
number of degrees of freedom is to assume that the constitutive coefficients are known 
throughout the domain, except at some locations where unknown inclusions may be 
present. The asymptotic expansions in the size of the inclusion mentioned above thus 
provide a very efficient tool to understand what may or may not be reconstructed from 
data with a given level of noise. 

For elliptic equations, the existing works on the subject, see e.g. [21 0], typically 
assume the parameters jumps across the interface of the inclusion. One of the main 
objectives of this paper is to consider the case of more general inclusions whose coefficient 
may vary at the small scale e and need not "jump" from the values of the background 
parameters. We also want to stress the similarities and differences between expansions 
for the diffusion equation V ■ D'^Vu^ = and the Helmholtz equation (—A + g^)^^ = 
with g*^ of order e~'^~^'^ for rj G [0,2]. In both cases of the diffusion equation and the 
Helmholtz equation when 77 = 0, we need to introduce local correctors and obtain a 
limiting influence at the domain's boundary that is non-linear in the parameters inside 
the inclusion. 

The maximal leading term in the expansion is always of order the volume of 

the inclusion. We construct expansions up to the order e'^'^. Going beyond this order 
of accuracy requires a more careful analysis of the decay properties of local correctors 
at inflnity than is available here, or the use of single and double layer potentials as in 
[2] in the case of constant coefficients inside and outside of the inclusion. Note that the 
cross-talk between two inclusions of volume 0{e'^) is also a term of order e'^'^. It seems 
therefore natural to stop the expansion at the order 0{e'^'^) for the influence of any given 
well-separated inclusions. 

Because our inclusions are modeled by somewhat arbitrary parameters that need not 
jump from the local value of the background parameter or are not constant, the limiting 
polarization tensors need not satisfy any property of positivity or deflniteness. On the 
contrary, we show that the polarization tensors vanish to flrst order for some types of 
inclusions, whose influence at the domain's boundary is therefore at most of order e'^^^ 
rather than e*^. Although we do not explore this aspect here, the proposed asymptotic 
expansions may be used to construct inclusions whose influence on the measurements is 
minimized in a prescribed manner. 

The rest of the paper is structured as follows. Section [2] is devoted to the derivation 
of the asymptotic expansions for the diffusion equation. The main tool in the expansion 
is a decomposition of the corresponding Green's function given in proposition 12. 1[ The 
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expansion obtained for smooth inclusions is presented in theorem 12.21 while the gen- 
eralization to more singular inclusions with possible discontinuities of the coefficients 
across the inclusion's boundary is given in theorem 12. 6[ We compare our expansions 
with those obtained in [2] for constant coefficients inside and outside of the inclusions 
in proposition 12.81 Section [2731 presents some properties of the polarization tensors that 
appear in the asymptotic expansions. In particular, proposition 12.121 shows that the 
leading polarization tensor vanishes for some non-vanishing diffusion coefficients inside 
the inclusion. Some proofs of the results are postponed until section HI 

Section [3] addresses local variations of the potential in a Helmholtz equation. The 
appropriate decomposition of the Green's function is shown in proposition 13.11 and the 
main result in theorem 13.31 The relationship between the expansions for diffusion and 
Helmholtz equations in regards of the Liouville change of variables is explored in section 
13.21 We show that the expansions in both settings agree up to order e'^"'"^. Most proofs 
are postponed to section HI 

2 Perturbations of the diffusion problem 

In this section, we are interested in the analysis of small inclusions in the diffusion or 
conductivity problem. As we have mentioned in the introduction, the reconstruction of 
diffusion or conductivity coefficients from boundary measurements is a severely ill-posed 
problem. One possible way to overcome this difficulty is to assume that the background 
diffusion coefficient is known and that the unknown part of the coefficient is localized 
and has small volume. 

Under such hypotheses, asymptotic expansions of the perturbed field in the volume 
of the inclusion have been derived in [B] when the inclusion is perfectly reflecting or 
insulating. These formulas have then been extended to more general inclusions in [5], 
and to higher orders in the volume and to domain with Lipschitz boundaries in [2]. In 
those references, the inclusion is modeled by a jump in the diffusion coefficient so that 
its first order effect on the boundary measurements is proportional to the inclusion's 
volume. The so-called polarization tensor contains the information about the inclusion 
that is available at this level of the asymptotic expansion. 

Such a setting for the diffusion coefficient prevents us from using the well-known 
change of variable q := that allows us to relate the diffusion equation to the 

Helmholtz or Schrodinger equation. Since one of the objective of the paper is to show 
the equivalence of the asymptotic expansions within the diffusion and Helmholtz frame- 
works, we first consider a regular inclusion without jump and derive the corresponding 
asymptotic expansions in section 12. 1[ We next generalize these formulas to the case 
with jumps in the diffusion coefficient in section 12.21 We also recover the formulas in 
[2] in the special case of constant coefficients in the background and the inclusion. Fi- 
nally, we present in section 12.31 some properties the polarization tensors involved in the 
asymptotic formula. 
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2.1 The case of smooth inclusions 



We consider the following system of equations: 

V ■ D^Vu^ = 0, in n, 

D'— = g, on dn, / u'da = 0, ^ ^ 

where f2 is a bounded open domain of dimension d >2 with Lipschitz boundary, a is the 
surface measure on dQ, and g e L'^{dQ) such that the following compatibility condition 
holds Jq^ gda = 0. It is assumed that D'^ is bounded from below by a positive constant 
independent of e and that -D^ satisfies the decomposition -D^(x) = -Do(x) + Di ( ) , 
where < Cq < Dq e C°°(n), Di G and vanishing in M'^\5, 5 being a 

bounded set with Lipschitz boundary. The properties of D"^ are summarized below: 

' D%x) > Co > 0, n a.e., 



Do(x), xGl]\xo + £:5, 

Do(x) + A(2^), xGxo + eB, 



We assume in addition that the domain of the inclusion is located away from the bound- 
ary in the sense that there exists do > independent of e such that 

dist{dQ,xo + eB) > do. (3) 

The Lax-Milgram lemma applied to ([I])-® yields a unique variational solution G 
H^{Q). Let us denote by U the solution with background diffusion coefficient Dq: 

V ■ DoVU = 0, in n, 

dU f (A) 

Do— = g, on dQ, / f/(x)rfa(x) = 0, ^ ' 

and introduce the related Green function G V{Q x Q) satisfying, for all fixed y in Q, 
Vx-/^o(x)VxiV(x,y) = -5(x-y), in fi, 

^ 9i^x^^_ 1 ^ on 91], [ iV(x,y)rfa(x) = 0. 

For all X G 12, the Lax-Milgram lemma yields again a unique variational solution U G 
and standard elliptic regularity results [7j implies that U G C°°(f2) since ^ 
C°°(12). We denote by F the fundamental solution of the Laplacian, namely 

-77-log|x|, d = 2, 

r(x) = <; 2^ 1 1 (6) 

d>3, 



{d-2)\S, 



\d-2 ■ 



where |S'(i_i| is the measure of the {d — l)-dimensional unit sphere. Throughout the 
paper, we use the following multi-index notations: for i = {ii, ■ ■ ■ ,id) & N"^, we define 
|i| = + . . . + i^, d'f = 9i7 ■ ■ ■ 9*7 and x^ = x\' ■■■x';}. We also define i\ = i^l ■ ■ ■ ij. 

One of the main tools in our asymptotic expansions is the following decomposition 
of the Green function A^: 
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Proposition 2.1 The Green function N can be decomposed, for (x, y) & Q x Q, as: 

iV(x, y) = /^o ^(x)r(x - y) + i?i(x, y) + R,{^, y) + R^iy), (7) 

where R3 G C°^(fi); for all y fixed in ^l, Ri{-,y) G W^'P{^1), with 1 < p < when 
d >3 and p < 00 when d = 2; and i?2(-, y) G H^{Q). Moreover, Ri is when x 7^ y, 
i?2 £ X and we have by construction that: 

VxiV(x, y) = Z}o^i(x) V.r(x - y) + Vxi?2(x, y). (8) 

ylfeo, admits the following asymptotic expansion for^ E B, y a.e. in dQ: 

VxiV(xo + £x,y) = 5^^Vx^9;iV(xo,y) + 0(£'^+^), (9) 

where 0{e'^~^^) denotes a term bounded in L'^{dQ) by Ce'^'^^ uniformly in x. 
Proof Let Ri he (uniquely) defined by 

Vxi?i(x,y) = ^|^r(x-y), / i?i(x,y)da(x)=0, 
^ok'^) Jan 



and i?3 be defined as 



|an|i?3(y) = - / '(x)r(x - y)da(x) 
Jan 



Since Dq > 0, Dq E C°°{Q) and F G L^qj,(M'^) for tlie values of p in the proposition, it 
follows that -Ri(-, y) G W^'''\Q). Moreover, Ri is C°° as soon as x 7^ y. In the same way, 
i?3 G C~(fi) since r(x) G C°^(M'^\{0}). We then verify that dZj) leads to ([HD and that 
plugging ([7]) into (jS]) leads to the system, for y E fl: 

Vx-I^o(x)Vxi?2(x,y) = 0, infi, 
c}i?2(x,y) 1 gr(x-y) /■ r ^ n 

which admits a unique weak solution thanks to the Lax-Milgram lemma since we verify 
that Jg^ij^^ + ^)'^cr(x) = 0. Since d^-§^i- — y) ^ L'^{dVl) for any multi-index (3 and 
y G f2, we deduce that d^R2{-, y) G so that elliptic regularity yields d^R2{-, y) G 

C°°(fi), and finally R2 G C°°(fi x 

Moreover, d^R2{-,y) is bounded in uniformly in y when y G fi' CC fi. To 

prove Q), we first remark from ([7]) that the trace c}^iV(z, y) [^^ is defined in L'^{dVt) 
uniformly in y when y E Vt' since i?i G C°°(f2\f2' x (9^i?2(-,y) G H^iyt) uniformly 
in y G ^2', and i?3 G C°°(f2'). This allows us to apply Green's theorem and obtain, for 
any (z,y) G fi' x 1], that: 



As y goes to dQ, the boundary integral converges for Lipschitz domains Q, see [2], to 

-J^ [ iV(x,z)da(x)-p.v / ^£[^iV(x,z)rfa(x) + iiV(y,z), (z, y) G fi' x 91^, 

where p. v. stands for the Cauchy principal value and the above integral operator is 
bounded in L?'{dVt). The first term belongs to C^iVt') and the second and the third terms 
to C°°(fi') with values in L'^{dVL). Using ([8]), this allows us to expand VxA^(xo + ex, y) 
and obtain ([9]). □ 

We first consider the case of a smooth inclusion by adding the hypothesis that 
Di is regular and compactly supported in that is Di G W^'°°{VL), with support 
suppi^i C B. In such a context, the trace of Di vanishes on dB. We have the following 
result: 

Theorem 2.2 Assume that Di G W^''^{VL) with support suppDi C B. Then the solu- 
tion to (CP-I^ verifies the following asymptotic expansion, a.e. on dVL: 



d d 



d_2+|i| + |j| 



\i\=i\j\=i 

d d d d d-2+\i\ + \j\ + \k\+l 

-EEE E Mg,,ffl^(xo) {d'D,^) (xo) g;iV(xo,y)|,, 

|i|=i |j|=i |fc|=oz=o,/+|fc|>o 

where M and are generalized polarization tensors given by 

= /z}i(x)V(x^+0°o(x))-Vx*dx, 2,JGN^ 
Jb 

Mf.^, = f (x) V(/.;.,(x) ■ VxMx, z, J, A; G / G N 

and t/ie functions are i/ie wmgne solutions in H]^^{W^) n C°°(M'^\S) to; 

' V-(Z^o(xo) + /^i(x))V05fc = -5°V- (Di(x)x'=Vx^) 

-«o(x„)E^^V.(«.,x)x™V^-H,,,), 



(10) 



m!(/ — Iml)! 
H=i ^ I 1^ 

[ 0;.,(x) = 0(|x|i-'^) 



as X oo. 



Here, 6i is the Kronecker symbol and the notation 0{e'^'^) in the expansion represents a 
term bounded in L^{dQ) by a constant depending on \\Di\\Loa and on \\g\\L'^{dn)- 

Remark 2.3 The function 0°^ solves the following equation in M*^: 

V-(Do(xo) + /^i(x))V<^;, = -V- (Z)i(x)x*^Vx^), 
0;.(x) = 0{\^\'~') 

so that 0^;, is computed from 0^, < m < I, iteratively. 
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Remark 2.4 We may recast the expansion in theorem 12.21 as 



d d d-2+|i| + |j| 

\i\=i\j\=i 

where the e-dependent tensor is given by: 

M^. = /" (x) V(x^' + ^j"(x)) ■ Vx'dx, z,jGN^ 

J B 

and the functions are the unique solutions in Hl^^{W^) fl C°°{W^\B) to 

V ■ (1 + Z}i(x)Do-i(xo + ex)) V^^^ = -V ■ {Di{^)D^\^q + £x)Vx^), 

^^(x) = C(|x|i-^) as |x| ^ cx). 

The asymptotic expansion of the theorem is then recovered by expanding \E'^ and 
Dq ^(xo + ex) in powers of e. 

There is another equivalent expansion to that of theorem 12.21 up to the order e^'^. 
We sketch its derivation in the case where Dq is constant. The right hand side of 
the equation for is equal to —D^^VDi ■ Vx-' — D^^Di^yJ . It turns out that an 
appropriate linear combination of Ax-' is of order 5'^'^^, so that we can replace in the 
definition of Mf- by $j solution to 

V ■ (1 + Di(x)Do-i) = -Do-VDi(x) ■ Vx^ 

$j(x) = 0{\^\^-'^) as |x| ^ 00. 

The appropriate linear combination is deduced from Af/ (xq + ex) = and from Taylor 
expanding U so as to obtain: 

d |j| 

= Af/(xo + ex) = A 5^ ^ 9^f/(xo)x^' + 0{e''+^). 

Remark 2.5 The leading order in the expansion is given by 

Yl M,,9W(xo,y)9^f/(xo). 
Kl=|j|=i 

The polarization tensor contributes only to higher orders. The polarization tensor 
M captures the correction when the background diffusion coefficient Do is constant in 
Xq + eB, whereas is the correction that needs to be added when Dq is not constant 
in xo + e-B. When Do is constant in xo + £-B, then M^^.^j = Mfj/^j^S^ so that the expansion 
then reduces to the classical formula: 

d d ii^2+\i\ + \j\ 

u%y) = ^(y) + EE ^ •!■! M,jd:,N{^o,y)d^U{^o) + 0{e"'). 

In this case, using the notation of remark 12. 4[ no longer depends on e and may be 
identified with 0^q. Note that the latter formula also holds when Do is non-constant 
away from the support of the inclusion xq + eB as remark 12.41 makes clear since only 
the values of Dq^ on the support of Di are involved. 
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The proof of the theorem is given in section |H Its main ingredients are the integral 
formulation of ([1]) and the decomposition of the Green function given in proposition [2]TJ 
Additional boundary effects, which are not considered here, appear at the order 0{e'^'^) 
when the geometry-dependent corrector i?2(x, xq + ey) of proposition 12.11 is expanded 
in powers of e. When Dq is constant, a proper factorization based on the technique of 
double layer potentials allow us to obtain arbitrarily accurate expansions; see [2]. 



2.2 The case of singular inclusions 



In the preceding section, we assumed that the perturbed diffusion coefficient was regular. 
We may generalize the above theorem to include the case where Di is in L°°(M'^) with 
support in B and with a possibly non-vanishing trace (if it is defined) at the interior 
boundary dB. This generalization is achieved by regularizing the singular perturbation 
so that we can use the preceding result and then by computing the limiting polarization 
tensors. We have the following result: 

Theorem 2.6 Assume Di verifies ([2]) with no further assumption on its interior trace 
on dB. Then admits the same expansion as in theorem \2.2\ with polarization tensors 
still given by f fTOj) . where 0^;, is now the unique solution in Hl^^{M.'^) fl C°°{M.'^\B) to: 



x G 



with 



as 



CO, 



V ■ (Z^o(xo) + /^i(x))V05, = -5? V ■ (Di (x) x^Vx^) 

,«-|m| / 



^-^ mUl — m ! V 

\m\=l ^ ' ^' 



X e 5, 



(12) 



^o(xo) 



9n 



^-^ m\[l — \m\v- 



\m\=l 



dn 
x) X 



l-\m\ 

m ' • 'jk 



dn 



Di (x) x^n ■ Vx^' 
X G dB, 



(resp. 



) denotes the 



Here, n is the outer normal to the boundary of B, 

outer (resp. inner) trace of on dB as functions in H^'^^dB). 

The proof of the theorem is postponed to section HI 

Theorem 12.61 has been proved in [2] by using single and double layer potential tech- 
niques when the background diffusion coefficient Dq is constant on the entire domain 
Vt and when Di is constant on B. Our result generalizes that of [2| to the case of non- 
constant Dq and Di for which layers techniques are not available. The first order of the 
expansion can also be obtained from the general formula proved in |1] and in [S]. 

Remark 2.7 The expansion in remark 12.41 still holds for singular inclusions with 
now the unique solution in H^^^iW^) n C°^(M'^\fi) to 







X G 



75, 



V ■ (1 + Di(x)Do"'(xo + ex))V^) 



-V- (A(x)Do~^(xo + £x)Vx^), 
C(|x|i-7 as |x| oo, 



X G 5, 
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equipped with the jump condition on dB: 



dm 

(l + Di(x)Z^o"'(xo + £x))^ =Di(x)Z}o^i(xo + £x)n- Vx^ ^edB. 



As in the end of remark 12.41 we could also derive a modified asymptotic expansion in 
the case of singular inclusions. 

The above asymptotic expansions are compatible with the slightly different expressions 
for the generalized polarization tensors obtained in [2]. We have the following proposi- 
tion: 

Proposition 2.8 Assume that Di is a non vanishing constant on B and that Dq is 
constant on the set xq + eB. Then u'^ verifies the following expansion, a.e. on dfl, 



d d 



an 



^ ^ pd~2+\i\ + \j\ 

Uiy)\an -EE ,,,, M,^'^(xo) 9iiV(xo,y)|,„ + 0{e''^ 



\i\=i\j\=i 

where Ai is the generalized polarization tensor given in JMl by 

Mij = Di [ n- V(x^+0j(x))x*ci(T(x), i,jeN'^. 

JdB 

The functions (pj are the unique solutions in Hlg^[M.'^) nC°°((M''/i?) UB) to the problem: 



Dn 



dn 



0, X G (RV^) U B, 



{Do + Di) ^ 
an 



cj>jiy) - ny)Do'D, / n ■ Wda{^) = 0{\y 

JdB 

The proof of the proposition is also postponed to section HI 



Din-Vx^ ^edB, 

^~'^), when |y| oo. 



2.3 Properties of the polarization tensor M 

In this section, we give some symmetry properties and estimates satisfied by the tensors 
M in theorems 12.21 and 12.61 

Proposition 2.9 Let Oi, I3i G M, where i belongs to a set a of multi-index I. Then, the 
polarization tensor M verifies the following properties: 



I Do(xo)A(x) I 



c?x < ^ a,a,Mi,- < I Di(x) | V ( a.x') 



(ix. 
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Proof. Using the definition of M, we fiave, 



and tlie system solved by 0% deduced from (fT2l) imply that, 



/ Di(x)Vx^ - V</.°orfx. (13) 
Jb 



Consequently, 

Y,c^^PjMi^= I A(x)v(5^/5,x^) • V(5^a,x^)dx 



Concerning item (ii), we remark from the above equality that: 

2 

a^x-^ ) j 

For the other inequality, we split the sum as: 



c/x. 



^aiajMij = / Di|v(^^a,x^)) ^rfx + / A V ^ aj0°o) ' v( ^ a^x^^rfx. 
Since Do(xo) +-Di(x) is strictly positive a.e. in Q, the Cauchy-Schwarz inequality yields 



[ Z}iV(5^a,0?o)-v(5^a.x^)rfx 



< ( / (Do(xo) + Z^i)|v(5^a,0°o 

i6/ 



(ix 



B Z)o(xo) + /^i 



a,x 



dx 



In the same way, equation ( fT3l) gives: 



' (Do(xo) + Z^i)|v(5^a,0;o 



J6/ 



rfx < 



B Do(xo) + 



c/x 



so that 



i,iG/ iG/ ^« uv uy I i 



Do(xo)A(x 



(ix, 



B Do(xo) + Di(x) 



a,x 



iG/ 



(ix. 
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This ends the proof. □ 



Item (a) of the proposition is very similar to the estimates obtained at the first order 
in [1] . Such estimates can be applied to verify the definiteness or not of the polarization 
tensor. In particular, it gives: 

I |2 f ^o(xo)£'i(x) ^-^ ^ \ \2 f n f \^ 

i5i;o(xo) + A(x) Jb 

so that for Di constant, M is positive definite when Di > Q and negative definite when 
Di < 0, as it was shown in [21 [5]. The only possibility to cancel the above sum is then to 
set Di = 0, which means that there is no inclusion. Therefore, an inhomogeneity with 
constant diffusion coefficient always generates a perturbation of order e"^ on the measure- 
ments. The situation is different when Di is not constant. Indeed, when Z}i(x)(ix < 0, 

then M is negative definite, and when r„ Tirr-^^TTTr^r'^x > 0, then M is positive defi- 

° ' J B Z)o(xo)+-Di(x) ' 

nite. But when Di{x)dx > while at the same time zw^^J+d7(x)'^-^ ^ ^^en M 
might not be definite for a suitable choice of Di as we now snow as an application of 
the intermediate value theorem. We show first that the functional Mij : — > M, 
Di Mij[Di] is continuous. 

Lemma 2.10 There exists a positive constant C, such that, for all finite multi-index i 
and j , we have: 

\M,,[D\] - M,,[Dl]\ < C\\D\ - DI\\l^(^b)- 

Proof. Take two perturbation D\ in with support in B and denote by 

M[D{], M[Z)^] the corresponding polarization tensors. Then: 

M,^[D\] - M,^[Dl] = I {Dl - Dl)V^^ ■ Vx^dx + / {Dl - Dl)V(j)%[Dl] ■ Vx^dx 

Jb Jb 

Introducing Wj := 0°o[-^i] ~ using the equations verified by [-^i] 

(p^QlDfl, we find the relation: 

/ (/^o(xo) + Dl)\\/Wj\^dx = - I {D\ - Dl)Vwj ■ (Vx^' + V (t)%[Dl])d^. 
Jud. Jb 

Since V0°o bounded in L'^{R'^), this yields the estimate 

||Vwj||L2(njd) < C\\Dl - L'iIIloo(b). 
Using ( fT4l) . we obtain the desired result. □ 

Lemma 2.11 There exists a perturbation Di G L°°{Q) with Di(x)(ix ^ 0, such that, 
for a given 1 < I < d, the component M(,^^e^[Dl] of the polarization tensor M vanishes, 
where ci is the l-th vector of the canonical basis ofW^. 
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Proof. Setting ctj = S^^ in item (ii) of proposition 12.91 leads to 



Now take a D\ such that JgZ}{(x)(ix < 0. Therefore, Mej^e; [-Dj] < 0. We then con- 
tinuously transform D\ into D\ such that DQ(^Q)+]y^) ^ ^ keeping D\dii. non 
zero in the transformation. Such a transformation exists: let indeed D\ be a bounded 
function in B with positive and negative parts -D^ and D\. We set D\(h(. > D\(hi. 
so that JgD{(x)(ix < 0. Letting the negative part D\ continuously go to zero then 
gives a possible transformation. For the resulting D^, we have Me;,e;[-Di] > 0. Since the 
functional Me,,eJ-Di] is continuous from L°°{B) to M, we deduce from the intermediate 
value theorem the existence of a Dl with j^D\dyi ^ such that Me^^gJ-Di] = 0. This 
ends the proof of the proposition. □ 
As a corollary of the previous result, we have 

Proposition 2.12 There exists a perturbation ^ Di & L°°{Q) with spherical symme- 
try such that Mij = 0. 

Proof. Consider an inclusion with spherical symmetry. We find that Mij = Mo6j 
when \i\ = \j\ = 1 so that the above lemma yields the existence of non- vanishing 
perturbation such that Mq = and consequently no term of order e'^ appears in the 
asymptotic expansion. □ 

The latter result is to be compared with the case where Di is constant for which there 
is always a contribution of order e"^ in the expansion provided the constant is not zero. 



3 Perturbations in the Helmholtz equation 

This section addresses the problem of small-volume inhomogeneities in the Helmholtz 
equation. As we did for the diffusion equation, we derive an asymptotic expansion of 
the perturbed solution in the volume of the inclusions. 

3.1 Asymptotic expansion and polarization tensors 

We consider the following Helmholtz (or Schrodinger) equation posed in a bounded 
Lipschitz domain Q of M.'^, d >2, and with < 5 for technical reasons: 

-Atr(x) + (go(x) + -l^^qJ^^^)) v%^) = 0, x G f], 

= ^ G L^idn) on dn, 

on 

where xq is a given point in Q, go ^ L°°{Q) is the background index or potential, 
and qi G L°°{Q) is a local perturbation, with support localized in a bounded Lipschitz 
domain B. We consider the case with only one inclusion, knowing that the results below 
generalize to the setting with several well-separated inclusions so long as the maximal 
order in the expansion is sufficiently small so that the inclusions do not interact at that 
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order. The perturbation has a magnitude of order e'^~'^, with rj G [0,2]. The most 
interesting case is 77 = 0, which corresponds to the strongest type of perturbation. The 
latter case allows to relate the asymptotic formula given in the preceding section to the 
one that we propose below for a particular form of the potential qi. 

When go is negative, the above system models waves propagating in a medium 
perturbed by a small inclusion of diameter e with a refractive index of order We 
refer to [9] and [8] for the case of high-frequency waves in dimension two perturbed by 
small inclusions with index of order one. The case go and gi constant with go negative and 
Tj = 2 has been treated in [2] with Dirichlet conditions instead of Neumann conditions at 
the domain's boundary. When go is positive, (fT5l) models e.g. diffusive light propagating 
in a medium with background absorption go and zones of different absorption coefficients 
in a small volume. The case rj = 2 has been investigated in dimension three in [3] for a 
constant background go and a constant perturbation gi. 

We denote by V the solution of the unperturbed equation 

-AV + qoV = 0, xGf], 

dV (16) 
— = 5( on dil. 
on 

When go = 0, we assume the normalizing and compatibility conditions: 

/ Vda = and / gda = 0, (17) 
Jan JdQ 

where a denotes the surface measure on dQ. According to f[T51) . this also implies: 

J Qi {—^) t'^(x)rfx = 0, when go = 0. (18) 

In order to obtain the existence and uniqueness of a variational solution to f|T6l) . we 
make the following classical assumption: 

(H-1) Let u G H^{n). Then 

/ Vn-Vi;(ix+ / qoUvdx = 0, for all G iJ^(fi), 
Jn Jn 

implies that u = 0. 

Under (H-1), an application of lemma S3] of the appendix yields a unique weak 
solution V G H^{Q) to f|T6l) . When go := 0, the same holds thanks to conditions f|T7|) . 
Since we need high-order Taylor expansions of V in the sequel, we make the additional 
assumption that the restriction of go to a neighborhood xo + eB' of the set xo + eB, 
with B CC B', belongs to C°°(xo -|- eB'). Using standard elliptic regularity [Tj and ([3]), 
we obtain that V G C°°(xo + eB'). When first order expansions are considered, then a 
L°°(f2) regularity for V is sufficient. Existence and uniqueness for (fT5|) uniformly in e 
for e small enough will be given in the sequel. When 77 g]0,2], no additional condition 
is required on gi. When 77 = 0, we add the following assumption: 
(H-2) —1 is not an eigenvalue of the bounded operator T defined as: 

T : L\B) ^ L\B), if ^ Tipiy) = [ gi(x)v.(x)r(x - y)rfx. 

Jb 
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Here, F is the fundamental solution of the Laplacian given in (H-2) is verified for 
instance when gi > a. e. in B or when the following RoUnick type [11] norm of qi is 
less than one, 

(v^^V^^^|r(x,y)|)'dxrfy < 1, 

for some p > 1, or when qi is a Bohm-like potential of the form 

Av/l + /^i(x) 

91 (x) 



for some C^(R'^) function Di with support in B such that 1 + Z^i > in 
The case d = 2 and 77 = is particular in the sense that 



so that we cannot expect the perturbation caused by the inclusion to be small in the 
general case. We thus need to add an additional hypothesis to be able to treat qi as a 
perturbation. It is the case under the following symmetry assumption: 
(H-3) When d = 2 and 77 = 0, we assume that the solution to f|T5l) verifies that 



j ^^° jt;-^(x)rfx = 0. 



Note that (H-3) is verified when e.g. qo = thanks to ( |T71) . We introduce the Green 
function iV(x, y) G V{Q x Q) of (fT6!) . which for each fixed y in Q, solves: 

-A,iV(x, y) + go(x)Ar(x, y) = 5(x - y), x G fi, 

giV(x,y) „ _ (19) 
— = on oil. 

When go = 0, iV has to be defined as in ([5]). is symmetric in its arguments. Hypothesis 
(H-1) is verified e.g. when qo > 0, Q a.e. (with the normalizing condition when go = 0), 
when go is constant and not an eigenvalue of the Laplacian equipped with homogeneous 
Neumann conditions, or when the following RoUnick-type norm of go is less than one. 




^AQ^\^AQM\\N{^,y)\Yd^dy < 1 



n Jn 



for some p > 1. We have the following proposition, which allows us to decompose as 
the sum of the whole space Green function F and a regular function: 

Proposition 3.1 We have A^(x,y) := r(x - y) + -R(x,y), where R{-,y) G H^{Q) n 
W'^'^{il') with p < when 3 < d < 5 and p < 00 when d = 2 for any f2' CC ^2 
uniformly in y & Q' . When go = 0, then R belongs to C°°(f2 x Q). Moreover, N admits 
the following asymptotic expansion for x G -B, y a.e. in dQ: 

VxiV(xo + £x, y) = J2 ^Vx^9;iV(xo, y) + Oie"^'), (20) 
where 0{e'^~^^) denotes a term bounded in L'^{dQ) by Ce'^'^^ , uniformly in x. 
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Proof. We consider only the case Q'o 7^ since the case qo = follows from proposition 
12.11 Plugging A^(x, y) := r(x — y) + -R(x, y) into f|T9l) leads for any y fixed in Q to the 
equation: 

- Axi?(x, y) + go(x)i?(x, y) = -go(x)r(x - y), x G 

dR{^,y) gr(x-y) (21) 
a = ^ ' 

anx 

Pick an y G fi' CC and for any v G H^{Q), consider the linear form: 

l{v):=- [ go(x)r(x-y)t;(x)dx- / ^^^"^ ~ ^\ ;(x)rfa(x). 

Then / is continuous in H^{Q). Indeed, on the one hand, r(x — y) is uniformly bounded 
for (x, y) G dQ x Q' which allows us to treat the second integral. On the other hand, 
r G Lfjjj,(M'^) with p < when d > 3 and p < 00 when d = 2 so that the Sobolev 
embedding H^{^1) ^ L'^{^), for q < when d >3 and q < 00 when d = 2 implies 

|/(tOI<c(||r|| + l)\\v\\m{n), 

for q' > when d > 3 and q' > 1 when d = 2, where Br is a ball of radius R large 
enough. Since > for < 6, we get the desired result. Note that for d > 7, 
the above linear form is not continuous as we may construct functions v G H^{Q) of 
the form |x|~" such that r(x)f(x) is not integrable in the vicinity of 0. Lemma 14.41 
then yields a unique -R(-,y) G H^[Q) uniformly bounded in y when y G fi' by choosing 
ao{u,v) = /q(Vm ■ Vv + uv)d'x and ai{u,v) = J^(go(x) — l)uvd:x.. Standard elliptic 
regularity [7] gives, for 1 < < when d >3 and p < 00 when d = 2, that: 

||i?(-,y)lkMn') < C{\\R{;y)\\HHn) + \\nLriBa)), 

so that R{-,y) G W^'P{n') uniformly in y G Q'. 

To prove fl20|) . we decompose R as i?(x, y) := i?i(x, y) + i?2(x, y) with 

-Axi?i(x,y) + go(x)i?i(x,y) = -go(x)r(x - y), x G fi, 

9i?i(x,y) (22) 
- = 0, on oiL 



dn 



X 



-Axi?2(x,y) + go(x)i?2(x,y) = 0, x G fi, 

gi^2(x,y) gr(x-y) (23) 

7S = ' 

anx anx 

Consider first (l22l) for y G According to lemma 1131 Ri{-,y) belongs to H^{Q) 
and is uniformly bounded with respect to y. Let B' be a neighborhood of B such that 
B CC B'. Since r(- - y) G C°°(xo + eW) uniformly in y G dO^and go G C°°(xo + 
we obtain from elliptic regularity that Ri{-,y) G C°°(xo + eB) uniformly in y G dQ. 
Now, i?2 is treated almost exactly as the term R2 in proposition 12.11 so we highlight 
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the differences. According to the previous results on R, the trace A^(x, z)|g^ exists in 
L'^{dQ) uniformly for z G fi' CC Q. Thus we have the following integral equation: 

i?2(z, y) = - / ^^(^~y) iv(x, z)rfa(x), (z, y)en'x Q. 
Jan dn^ 

As y goes to dQ, the integral converges to 

_p.v / ^£^:iZliV(x, z)da(x) + iiV(y, z), 
Jdn on^ 2 

where p. v. stands for the Cauchy principal value and the above quantity makes sense 
in LP'[dVL) uniformly in z G ^2' so that i?2(z, ■) G L'^{dVL) for all z G Moreover, we 
verify that i?2(z,y) satisfies in the distributional sense, for z G y G dVL, 

-A^i?2(z, y) + go(z)^2(z, y) = 0, 

so that we conclude from elliptic regularity that _R2(-,y) G C°°(xo + sB) with values in 
LF'iydVl). Classical Taylor expansions then yield fl20|) . □ 
We come back to f|T5l) and state the following result. 

Proposition 3.2 Assume that (H-2) is satisfied when rj = and (H-3) is satisfied 
when d = 2 and t] = 0. Then, under assumption (H-1), there exists Sq > 0, such that 
for all < e < Eq, the system l[T5\) admits a unique variational solution v'^ G H^{Vi). 
Moreover, the restriction of f ^ to the set Xq + sB verifies the following decomposition 

w"(xo + ey) = y(xo + ey) + e'^^'iy) + 6'^-'^+" r%y) + 0{e'^'^^), y a.e. in B, (24) 
where ^'^(y) := Xlj^-fio ^^"'^(^o)0j(y) (^^'d (fij is the unique solution in H^{B) to 

<f)] + e^T<f)] = -T^^, yG5, (25) 
and the unique solution in H^{B), fory G B, to 

r'{y) + e'^Tr^y) = [ qi (x) v'{^o + ^x) (i?(xo + ex, xq + sy) - (27r)-^ loge) rfx. 
Jb 

The operator T is defined in (H-2) and the function R in proposition \3.1\ whereas 6^ 
is the Kronecker symbol. The notation 0{e'^'^'^) represents a term bounded in H^{B) by 
Ce'^'^^ . The remainder is bounded in L'^{B) independently of e when d = 3, by Ce~"^ , 
for any a > when d = 4, by Ce~^ when d = 5, and by C \ loge\ when d = 2. When 
d = 2 and rj = 0, then is of order 0{e) thanks to (H-3). When q^ = 0, then is 
bounded in H^{B) independently of e for any d. 

We then have the following theorem: 

Theorem 3.3 Under the hypotheses of proposition ISTB, the solution to to ( [73]) sat- 
isfies the following asymptotic expansion, almost everywhere on dQ: 



d+l d+l d~2+r,+ \i\ + \j\ 

v^iy)\an = ny)U-EE -vfi {Q^, + e^Ql)^^V{^o)^^N{^o,y) 

\j\=o\i\=o 



\an 



-£2('^-2+'')/^(y) + C(£ 



2d\ 
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where 0{e'^^) is a term bounded in L'^{dVL) by Ce^"^ and for {i,j) G N*^ x N'^, 
Q,, = I gi(x)x^xMx, Ql = I gi(x)</.^^(x)xMx, 

J B J B 

riy) = / gi(x)r^(x)A^(xo + 5x,y)dx. 
Jb 

The remainder is of order: 0{ \ logel) when d = 2; 0{1) when d = 3; 0{e~°') 

for any a > when d = 4; and 0{s~^) when d = 5. 

The proofs of the proposition and the theorem are given in section When 77 > 0, 0J 
still depends on e. We may then expand the operator {I + e^T)~^ in terms of Neumann 
series up to the right order. We include the term in the formula because we need 
its explicit expression below to make the link between the asymptotic expansion for the 
diffusion equation and that for the Helmholtz equation. 

In the particular case where go constant and positive, 77 = 2, gi is constant, and the 
inclusion is centered at Xq so that xc/x = 0, we find for d = 3 that 

^^(y) = V{y)-e'q, (^J^{l + e'^l)d^^V{^o)N{^o.y) 

- ^1 E E I"! ( / ^'^'^^) ^'^(^0' y)^'^(^o) + 0{e'). 

According to 0q verifies 0q = — Tl + 0{e^^ so that we recover the asymptotic 

expansion given in [3]. 

The tensor Q is clearly symmetric. When q\ is constant and not identically zero, there 
is always a contribution of order -^^ expansion, while for spatially varying gi, 

the first order contribution can vanish for instance by choosing gi such that girfx = 0. 

3.2 Relation between the diffusion and Helmholtz equations 

We now compare the asymptotic expansions for the solution to the diffusion equation 
([1]) given in theorem 12.21 and for the solution to the Helmholtz equation f|T5l) given in 
theorem 13. 3[ It is well-known that a solution to the diffusion equation 

V ■ D^u = 0, 

with D G C^(R'^) for instance and strictly positive, also satisfies a Helmholtz or Schrodinger 
equation of the form 

Our purpose here is to verify that the polarization tensors obtained in the diffusion 
and Helmholtz frameworks are indeed the same for the specific form of the potential 
gi that allows us to transform one equation into the other. As in section [2l we define 
D'^ix.) = -Do(x) + Di{ ^~^" ) and to simplify the presentation, assume that Dq is constant 
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in Q. We assume that Di G C^^Q) with support included in B and that Dq + Di is 
strictly positive in Q, so that we can define 



Ay/Do + Dl(x) 

a/ Do + A(x) 

We then consider the function which satisfies f|T5l) with go = 0, = and gi defined 
as above. With such a choice, the quantity 



Do + Di 



• X-XQ ' 

s 



solves ([T]). Since 77 = 0, we may expect from the expansion given in theorem 13.31 that 
the inclusion induces a correction of order e'^"^ whereas the same inclusion induces 
a correction of order e"^ in the diffusion equation. Some simplifications due to the 
particular form of the potential qi must render the correction of order e'^ in the Helmholtz 
framework as well. We state the main result of this section: 



Proposition 3.4 When qi has the form ^EMj, then we have 

d+l \j\ 

EV^'^(^o)(go. + gy = o{e'^% (27) 



d+l |j| 

^^9W(xo,y)(g.o + Q°o) = 0{e'^'). (28) 



i=0 



Here, the index of the polarization tensors represents the vector ofN'^ with components 
all equal to zero. We have the following relation between the polarization tensor M in 
the context of theorem \2.2\ and the polarization tensor M := \/Dq{Q + Q^) in the context 
of the Helmholtz equation: 

Mi^=M,^, |^| = |j| = l, (29) 



iji=iiii=i 



The proof of the proposition is given in section 14.21 Equations (127|) and (128|) imply 
that the two first orders in the expansion of theorem 13.31 vanish so that the correction 
is of order e*^. Equations (l29l) and (l30l) show the equivalence of the tensors Mij and 
Mij for 1^1, 1 j| < d + 1 up to an error of order e'^"'"^, which is sufficient to show that the 
asymptotic expansions on and agree up to the order e"^"^. The proofs can in fact 
be modified to show the equivalence at higher orders as well, e.g., for any r G N, 

^ pKI+lil 
\j\=i\i\=i 

Furthermore, denoting by {rriij) the modified polarization tensor obtained from $j at 
the end of remark 12. 4[ we can show in this context the strict equality between the 
Helmholtz and diffusion tensors, that is Mij = rriij, for all i,j. 
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4 Proofs of the main results 



4.1 Asymptotic expansions for the diffusion equation 

We now prove theorems 12.21 and 12.61 and proposition 12.81 

Proof of Theorem \2.S[ The starting point of the proof is the formulation of ([T]) as 
the following integral equation: 

^^(y) = U{y) - [ D, f"^^) Vw^(x) ■ VxiV(x,y)c?x, 

JxQ+eB V ^ / 

= f/(y) -e'' [ (x) Vn^(xo + ex) ■ VxiV(xo + ex, y)dx. (31) 
Jb 

The above equation is justified rigorously as in the derivation of (1751) in lemma I4l2l of 
the appendix. We highlight the main differences. According to proposition 12. 11 we have 
VxiV(x,y) = Do-^(x)Vr(x - y) + V,i?2(x,y), with W^R^i^y) G L\n) for every y in 
Q so that the above equation makes sense in L'^{Q) and therefore almost everywhere in 
Q thanks to the Young inequality since Vm^ G L^(f2) and VF G Ll^^{M.'^). The integral 
equation fl3Tl) is obtained from the variational formulations of ([1]) and (jlj): 

/ D'Vu'-Vvdx = I gvda{-K) = ! D^VU-Vvd^, (32) 
in Jan Jn 

for all V G H^{n). Then, let (p G ^^(fi) and set w(x) := J^N{x,y)(p{y)dy. Thus v is 
the unique solution in H^{Q) to —V ■ DqVv = ip equipped with homogeneous Neumann 
conditions and the normalization jQ^vda{x.)=0. As in the proof of fl75|) or in the proof 
of proposition 12.81 we verify that Fubini's theorem applies and that 

^ (^^Do(x)V^i(x) ■ VxiV(x,y)dx-M(y)^ ^(y)rfy = 0, Wu e H\n). 

Applying the latter equality to both and U, gives (l3T|l together with (132|) . 

To continue the proof of theorem, we write u"^ = U as the sum of the unperturbed 
solution U and a corrector w^, solution of 

V ■ (Z^o(x) + Di (^)) Vw' = -V ■ Di (^) Vf/, in fi, 

^ = 0, ondn, / M;"(x)c/a(x) = 0. ^ ' 

Since both and f/ belong to H^{VL), then w;'^ G iJ^(r2) and we deduce from (1551) that: 

for some xq + SqB d Bq dd VL with eo > so that, from standard elliptic regularity, 

||Vw;^(xo + e-)IU2(B) < C'||^i||L°°(B)||Vf/||L°o(Bo) <C\\Di\\L^^B)\\9\\L-2{dQ), (34) 

for some constant C > 0. We need an approximation of the corrector w"^ up to the order 
e'^ and so that we decompose it as w'^(xo + ex) = \E''^(x) +r^(x), where is a remainder 
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of order e'^ in a sense made precise below. Finding an asymptotic expression for then 
amounts to calculating \I^^(x) and showing that is indeed of order e'^. To this aim, 
we use ( 13T|) to obtain an integral equation for w'^ verified a.e. in Vt: 

w'{y) = -e^ I Di (x) V [w' + U] (xq + ex) ■ VxA^(xo + ex, y)rfx. (35) 
Jb 

We then decompose A^(x, y) following ([8]). Plugging ([8]) into fl35l) . setting y := xq + ey 
for y G -B, and using the homogeneity Vr(£:x) = £:^~'^Vr(x), we find 

w'{^o + ey) = -£ / Z^i(x)Z}o^i(xo + £x)V[«;^ + f/](xo + £x) ■ Vxr(x-y)dx 
Jb 

-e'^ I Di (x) V [w'' + U] (xo + ex.) ■ Vxi?2(xo + ex, xq + ey)dx. 
Jb 

We shall prove that the contribution involving R2 above is of order 0(e'^) and that up to 
an error of the same order, we may replace Dp ^(xo+ex) and [/(xo + ex) by Dq^(xo+£:x) 
and Ud{xo + ex), respectively, where for H = Dq^ and H = U, we have defined the 
Taylor expansion to order d: 

d \m\ 

H,{xo + ex) = V — (d'^H) (xo) x'". (36) 
^ — ' ml 

\m\=0 

Note that eWw^xo + ey) = V^'^(y) + Vr^(y). We thus want ^'^(y) to solve: 

^^(y)+To,rf^^(y) = -£ /z}i(x)Z}o^^(xo + £x)Vf/d(xo + £x)-Vxr(x-y)cix, (37) 

Jb 

where we have introduced the notation 

To,rf^(y)= / Z}i(x)D-i(xo + £x)V*(x)- Vxr(x-y)rfx. (38) 
Jb 

The above equation is the integral formulation of 
A^- + V ■ {D,{x)D,l{xo + £x)) V^^ = -sV ■ (Di(x)Do;^(xo + £x)(Vf/rf)(xo + ex)). 
We now thus expand ^(xq + ex) in the definition of Tq^^i to obtain: 




To^(y) := / A(x)Do"'(xo)V^(x)-Vxr(x-y)rfx. 
Jb 

Expanding Ud and Z^q^ in (!37|) . and setting 

d d |j| + |fc| 

^^{y) = Z^o(xo) E W (^'f^) (^0) (^'^0 ^) (xo) *,\.(y), 

|i|=l|fc|=0 
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leads to the following equation for 




-Doi^or' f Di(x)x'=Vx^ - Vxr(x-y)rfx. 
Jb 



Equating like powers of e, we verify that ^I/jfc(y) = X]f=o 7r*i^ifc(y)' where ^^-^i, solves the 
following integral equation a.e. in every bounded set of M'^: 

(/ + r„)^;,(y) = - t / C,(x)x-V^;^H(,).vj(x-y)rfx 

|m|=l ' 

-5?Do'(xo) / /^i(x)x'^Vx^ -Vxr(x-y)(ix. 
Jb 

Existence and uniqueness of solutions in Hl^^iW^) fl C°°(M.'^\B) to the above equations 
follows from lemma \^72\ of the appendix: we first prove the result for 0^^, then for 
which depends only on 0°^, and finally for all 0^ iteratively. Moreover, according to the 
lemma, 0^-^ solves the system of differential equations given in ( fTTI) . The function "if^ 
thus belongs to the space Hl^^{R'^) n C°°(M'^\S) by construction. When Dq is constant 
and equal to -Do(xo) in the set xq + eB, we do not need to expand Dq^. We thus have 
Do} = Do\^o) and Q can be identified with \&^q. 

We then verify that the remainder r^(y) = w^(xo + ey) — ^E'^(y) belongs to H^{B) 
by construction and moreover solves the integral equation: 

(/ + To,,y{y) = S'iey) - e''+^ I Di(x) [Si(x)Vw;^(xo + ex) + S2(x)] ■ Vxr(x - y)(ix 

Jb 

where Si is the remainder of the d + 1 order Taylor expansion of D^^ (xq + ex) (so that 
Do ^(xo+ex) = i:)^^(xo+£:x)+e'^+^S'i(x)), S2 the remainder of ^(^o+e^x) Vt/(xo+£x) 
and where we have defined 

S^{ey) = -e"^ I Di (x) Vm"(xo + ex) ■ Vxi?2(xo + ex, xq + ey)rfx. 
Jb 

We may now decompose r' as r"(y) := rf (y) + r|(y) + ^^(ey) with: 

(J + To>f(y) = -e /Di(x)D-i(xo + ex)V5^(ey)-Vxr(x-y)rfx, 

(J + To^dKiy) = -e"+2 /■ [Si(x) Vw7^(xo + ex) + S2(x)] ■ Vxr(x - y)rfx. 

Jb 

We know from the hypotheses in ([2]) that for all y E B, -Do(xo + ey) + -Di(x) > Cq > 0, 
so that setting < e < eo for eo small enough, we have 1 + -Di(x)Z}q ^(xg + ey) > Ci > 0, 
for another constant Ci independent of e. An application of lemma 142) then yields that 



21 



rf and r| are uniquely defined in H^^^iW'^) r\C°°(R'^\B). Moreover, following lemma W?2\ 
we have the estimates: 

II Vrf ||^2(iRd) < C£||Di||loo(b)||VS'^(£-)||l°°(b), 

||Vr^|U2(Kd) < C£'^+lDi|U^-(B) (II Vw^(xo + E-)\\l2(b) + ||I^o"'Vf/||crf+i(Bo)) , 

< Ce'^^'^\\Di\\1^^^^\\g\\L2(^QQ), 



according to (IMI) and by elliptic regularity, where -Bq is as above ( iMll . It thus remains 
to estimate S*^. From proposition 12. ![ we know that R2 G C°^(r2 x Q), which yields: 

||V5^(£-)l|L-(i?) < Ce^PlllL-(B) (||Vti;^(xo+£-)llL2(B) + ||Vf/|Uoo(Bo)) X 

X II VxVy-R2||L°=(BoX_Bo) < ^"'l|-Dl||ioo(5) ||5'||L2(an)- 

Gathering the different estimates for rf, and S**^, we obtain that 

To conclude the proof, we go back to ( 13T]) and take the trace on dQ. Plugging Vw^{'Xq + 
£x) = e"i(V^'^(x) + Vr^(x)) into ([31]), it just remains to expand Vf/(xo + ex) G C°°(5) 
and VxA^(xo + ex, y) thanks to ([9]) since we find, a.e. in 9^2, that: 

^%y)\on = Uiy)\on " ^'/^^i W (^^^(^o + ^x) + lvvl/^(x)) ■ iV(xo + ex, y)|,^ dx 

The asymptotic expansion of remark 12.41 is obtained by decomposing w'^ slightly differ- 
ently. We write tf^(xo + ex.) = ^'^(x) + r^(x), where \E'^ is now given by 

*^(y) + r^*^(y) = -s [ D, (x) D,\xo + £x)Vt/,(xo + ex) ■ Vxr(x - y)rfx, 

Jb 

T^*(y) := /di(x)Do-1(xo + £x)V^(x)- Vxr(x-y)dx. 
Jb 

We then verify that the remainder is of order e'^ and that expanding Ud and setting 
^^(x) = ^|^^.|^^ ^ (d^U) (xo) leads to the desired result. □ 

Proof of Theorem \2.6l Let Di be a non-regular perturbation and let be the 
cut-off function with support in B defined as 

r X^(x) = 1, for X G -B such that dist(x, (95) > rj, ^^^^ 
\ x''(x) = 0, otherwise. 

The parameter 77 will be adjusted according to e. Let now p''(x) := ri~'^p{ri~^x) be a 
standard moUifier and let := * (x^'Di). We verify that G C°°(]R'^) and that 
its support is included in B with a vanishing and continuous trace at the boundary. 
We can then apply theorem 12.21 to obtain an asymptotic expansion for the solution 
associated to D^. Since the error term of order e'^'^ depends only on ||D^||j;^oc.(-]jjd) - which 
is bounded by ||Z}i||^oo(-]Rd) - it suffices to look at the limit of the different polarization 
tensors to find the limiting asymptotic expansion. 
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Since -Do(x) + Di is bounded from below by Cq, this property is still verified 

by the regularized diffusion coefficient so that, according to fl7^ of lemma 14.21 of the 
appendix, the function 0^'^ associated to satisfy by induction the estimates, for 



for any bounded set A. This yields that V<p)k bounded in L'^(M'^) independently of r] 
and so is (p^J! in H\A). Defining the set E := {{j,k) G G N, < \k\,l < d} 

with cardinal \E\, we may thus see {^j'^js as bounded in {H^{A)y^^ and extract a 
subsequence as r/ ^ converging strongly in {L'^{A)y^^ and with gradient converging 
weakly in (^^(M'^))!^! to a limit {^Jfcjij. We obtain that G i/iclK'^) and V^J^ G 
L^(R'^). To find the equation solved by V^'^, we consider the weak formulation verified 
by 05.'^, which is, for all functions if G Hl^M'^) such that V(p G ^^(M^), ^ 
as -R — *■ 00, where Sr denotes the sphere of radius R, 

[ (Do(xo) +Z^?(x))V05'fc^- V¥5cix= / D'>{^)^^W -Vipd^ 
Jm.'' Jb 

|m|=l ' 

The above formulation is justified in lemma W7I\ below: see fl77l) . Since D\ converges 
strongly in any L^iW^) for 1 < p < cxd, we can pass to the limit in the non-linear terms 
above and obtain the following limiting equation: 



(40) 



/ (Z}o(xo) + /^i(x))V05fc- V<^dx= -5° ! Z^i (x) x'= Vx^' ■ Vv? rfx 
-^o(xo) E ^u/l\m\v i^'"^o') (xo) / A(x)x-V0;-l-l(x).V^dx. 

\m\=l ' ''^ 

To obtain the behavior of </)^-^ at infinity, we use the integral formulation given in flTSl) 
of lemma l42l for the subsequence 0^'^ and obtain, a.e. in every bounded set Q' C M°': 

(/ + To)0;,^(y) = - E / ^inx)x-v4H,.(,).vj(x-y)rfx 

|m|=l ' 

-SfD^\^o) [ D'1 (x) x^Vx^' ■ Vxr(x - y)dx. 
Jb 

The above equation makes sense in L^(fi') and therefore almost everywhere in Q' since 
V(^5'^ G / = 0, ■■■ ,ci, and Vr G Aocl^'^) so that the right hand side is fi- 

nite thanks to the Young inequality. Consider now a compact set K G such that 
dist(ii', B) > C > 0. The above equation is then verified uniformly in K and moreover 
e C^{K). Since V0^'^ converges weakly to V^^-^ for < / < (i and converges 
strongly, it follows from the above equation that 0''^ is a Cauchy sequence in C^{K) so 
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that it converges uniformly to the solution, for all x G -ft', to 



E 

|m|=] 



m\{l — \m\)\ 



Di (x)x' 



B 



x) ■ Vxr(x - y)(ix 



-S^D^\^o) / /^i(x)x'^Vx^ - Vxr(x-y)rfx. 



(41) 



B 



The fact that Vxr(x 



y) = Cdyl^^"*) for X e 5 and y E K yields that ^^^(y) = 
0(|y|^~'^) for such values of y. It is then not difficult to see that pUj) is the weak 
formulation of the problem given in the theorem. Notice that equation fHTl) is also valid 



a.e. in A since (pjf^ G L'^{A) for any bounded set A. Uniqueness follows from flKIl) 
and the behavior at infinity: the right hand side of (HOj) vanishes when we consider the 
difference of two possible solutions. Since those solutions are sufficiently regular, taking 
that difference as a test function implies the difference is a constant which must be equal 
to zero according to the vanishing limit at infinity. 

Now that we have the expression of the limiting 0^-^, it suffices to pass to the limit in 
the polarization tensors using the weak convergence of V(/)^'q and the strong convergence 
of and to choose rj small enough such that all the errors terms coming from the 
different passages to the limit are smaller than Ce'^'^. □ 



Proof of Proposition \2.S\ When Dq is constant on the set xq + eB, only the sum 
involving the polarization tensor M remains in theorem 12.21 as we have mentioned in 
remark 12. 5[ We thus start from the expression of M given in theorem 12.61 and define 
fj := 4>^Q — (j)j. A proof of the existence and uniqueness for (pj can be found in |2]. Owing 
the definitions of </)°o and 0^, we find that fj e HlM'^) n C°°(M'^/S) by construction 
and is the unique weak solution to 

V • (Z^o + lIi?(x)A)V/,- = -]I,5(x)Z}iAx^ xgM^ 

equipped with the condition at infinity: 



/,(y) + T{y)D,'D, [ n ■ Vx^rfa(x) = 0{\y\'- 

JdB 



Here, ll^ is the characteristic function of the set B. When |j| = 1, we obtain fj = 0. 
When Ij'I > 2, we need to sum over j to show that fj is small in an appropriate sense. To 
this aim, we derive an integral equation for fj from that of As we mentioned 

in the proof of theorem 12.61 (HTj) is verified a. e. by (t>^n so that we have 



DoD^'^liy) 



V0°o(x) ■ Vxr(x - y)dx - / Vx-'' ■ Vxr(x - y)dx. (42) 



B 



Since (pj is harmonic in BUM.'^\B, we deduce from elliptic regularity in Lipschitz domains 
(see e.g. [2]) that (pj G H2{B) so that its inner normal derivative at the boundary dB 
belongs to L'^{dB). This allows us to express (pj in terms of single layer potential, using 
the jump of its normal derivative at the boundary given in proposition 12. 8[ as 



Do^rV,(y) 



dB 



dn 



x) + n ■ Vx-'' r(x - y)da{x). 



(43) 
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The latter equation is verified in L^{A) for any bounded set A C M.'^, and tlius a.e. since 
||r(x — is uniformly bounded in x G dB. Moreover, since (f)j is harmonic in B, 

we have for any G H^{B): 



B JdB dn 



(f (icr(x). 



Let ip G and set = J^'?/;(y)r(x — y)(iy. Using the Young inequality and the 

fact that r and VF belong to Ll^^{M.'^), we verify that ip G H^{B) so that it can be used 
as a test function. Moreover, to be able to use the Fubini theorem, we apply as in the 
proof of lemma the Sobolev inequality 14.31 to conclude that V<;/'j(x) ■ VF(x — y)^/'(y) 

belongs to L^[B x A). In the same way, ^ (x)F(x — y)'>p{y) belongs to L^(dB x A) 



smce 



dB J A 



dn 



(x)F(x-y)^(y) 







dn 





da{x.)dy<C ^— l|J- llSall''/^llL°°(A), 

- L^{dB) 

for a ball of radius a large enough. We may thus write: 

F(x-y)da(x) ) ^/;(y)dy = 0. 



V0j(x) ■ VF(x-y)rfx- 

A \JB JdB 



dn 



Plugging (143!) into the latter equation yields: 



DoDrV,(y) + / V0,(x) ■ VF(x - y)t/x + / n ■ Vx^' F(x - y) rfa(x) ij{y)dy = 0. 
A \ Jb Job / 

Integrating fj42|) against ip, subtracting the equation above and performing an integration 
by parts, we find: 

^ (^DoD^'f,{y) + ^ V/,(x) ■ VF(x - y)dx - ^ Ax^' F(x - y) rfx^ ^{y)dy = 0. 

The quantity under parentheses belongs to L'^{A). Thus, we deduce by density that 
the above relation holds also for any ip G L'^{A) so that fj solves the following integral 
equation, a.e. in every bounded set Q' C M*^, 



DoDi'f.iy) 



B 



V/,(x) ■ VxF(x - y)dx + / Ax^F(x - y)rfa(x). (44) 



We now show that an appropriate linear combination of the //s is of order e'^'^^. First 
since Dq is constant in xq + eB, At/(xo + ex) = for x G i? according to (j4]), so that 
using the notation in ([36D, we get that At/d(x) = A([/d(x) - f/(xo + ex)) = 0{e'^'^^' 
uniformly in B. As a consequence, we have 



i?^(x) := A?7,(x) = J2 — 9^f/(xo)Ax^- = 0{e'+'). 
|j|=i 
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Thus, defining 

d 

e 



lil=i ■'' 



it follows: 



\3\ \3\ C 

lil=i lil=i "^^ 

= Z^i V ^d^Ui^o) [ V(x^' + 0j(x)) ■ Vx^dx + Di [ VF"(x) ■ Vx^c^x. 



lil=i 

According to the definition of G i/iocl^^^) n C°°(M'^/fi) solves: 

V-(Z^o(xo) + ]Ib(x)Di)VF^ = -]IB(x)Dli^^ xgM^ (45) 
equipped with the condition at infinity: 

F^(y) + r(y)Di / i?^dx = 0(|yr-'^). (46) 
Jb 

Following (jHj), thus solves the integral equation, a.e. in every bounded set Q' C M*^: 

F^(y) = -DiD,\^o) f [VF^(x)- Vxr(x-y)-/?^(x)r(x-y)]rfx, (47) 

Jb 

so that Young's inequality gives 

WF'WmB) < C\\VF'\\mB) + C\\R'\\l2^b) < C||VF^|U2(5) + 0(£<^+i). (48) 

Let Br be the ball of radius R with B GG Bji and denote by 5^ its boundary. Inte- 
grating 0451) on Bfi against leads to 

(Do(xo) + ]lBF'i)|VF"p(ix= -Di / i?"F"rfx+ / -—F'da{^), (49) 



Bu JB J5« 5n 

where a is the surface measure on Sr. We may recast condition (H6il using the integral 
equation (147!) for F"^ and its derivative as 

d-F'iy) + a"r(y)A / R'd^ = O {{\\VF^lhb) + ||i?1|L2(B))|yr^'-l°l) , (50) 

Jb 

for a multi- index a with |a| < 1. Consider first d > 3. Then VF^ G L^(]R'^) and the 
boundary integral in (H9!) goes to zero as i? tends to infinity so that 



B 



(Do(xo) + nBDi)\VF'\^d^ = -Di I R'F'd^. (51) 

i 

This yields, together with (HHj) : 

||VF^||i.(j,,) < 0(£2('^+^)) + ||i?^|U.(5)||VF^|U.( 
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so that II VF^||i2(Rd) = 0{e'^^^). Consider now the case d = 2. We cannot use the same 
approach since does not vanish at infinity. Using (l50l) for y G 5*^, we have 

|F^(y)| < c((logi? + i)||i?^|U.(s) + i||VF^|U.(B) 



dn 



(y) 



L^B) + II VF^ 11^2(5) 



SO that, since ||i?=||^2(B) = 0{e'^+^ 



I ^F^dai^) < Ce2('^+^)+C«.'^+iVF^|U2(^) + §||VF^||i2 
Jsr on ' ' R-^ 



Since, according to hypothesis [21 Dq^-xq) + I[_b-Di > Cq > a.e. in M*^, it follows from 
( l48l) . (I49l) and the above inequality that: 



CollVF^III^ < Ce^C^^^) + + r/)||VF^||i2 



for any 77 > 0. It suffices finally to set rj small enough and R large enough so that 
+ ?7 < Co to obtain 



||VF^|U2(5)< ||VF^|U2(B,) = 0(£ , 
We end the proof with the following integration by parts: 

lil=i ^' \j\=i ^' •'^^ 



-Di / i?^(x)xVx + Di / VF^(x) ■ Vx'dx, 
Jb Jb 

d „ 

Di V i-9^f/(xo) / n ■ V(x^' + 0,(x))x'rfa(x) + 0{e 

m=i -^^^ 



lil 



d+U 



|J|=1 



J! 



which shows that the error terms generated by M and A4 agree up to an order 0{e'^+^). 
□ 



4.2 Asymptotic expansions for the Helmholtz equation 



We now prove proposition 13.21 theorem 13.31 and proposition 13. 4[ 

Proof of proposition \3.2i We write v"^ := V + so that the corrector satisfies: 



dn 



-Aw^ix) + ( go(x) + ( 1 1 wHx) 



r2-r, 



e 



0, on dVL. 



(52) 
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We need to show the existence of w'^. We first show the existence and uniqueness of a 
solution to the integral formulation of fl52l) . which formally reads, y a.e. in VL: 



TV(y) = / gif^^^')v^(x)iV(x,y)dx. 

J^o+eB \ ^ J 

We consider first the case c/ > 3. Using the decomposition of given in proposition 13. II 
and denoting by w* the restriction of w^(xo + £:y) to B (we do not write the dependence 
on e to simplify), we recast the above system as 



w* + + e'^-^+'^R'w* = -T'Vi^o + ey), y E B, 

Tw*{y)= / gi (x)u;*(x)r(x-y)cix. 



(53) 



R^w*{y) = / gi (x) w*(x)_R(xo + ex, xq + ey)(ix. 
Jb 



We have used the homogeneity r(ex) = £:^~"'r(x) when d > 3. Since T and are 
compact operators in L'^{B), they have discrete spectra. Indeed, since F, VF G Lj'^j^(R'^), 
we have, using the Young inequality for any ip G L'^{B), 

\\T(p\\l2{B) < lklV'||L2(B)||F||/,i(B„) < |kl||L°°(B)||'/^||L2(B)||F||Li(B„), 

where Ba is a ball of radius a large enough. Thus, proceeding analogously for VT, 

\\T\\c{L'2(B)) < lkl||L°°(B)||F||Li(B„), II Vr||£(L2(5)) < ||gi||L°^(B)||VF||ii(B^), 

and compactness stems from the Rellich theorem. The same holds for since it is 
Hilbert-Schmidt as -R(xo + £:-,xo + e-) belongs to L^{B x B) (though not necessarily 
uniformly in e; see below) according to proposition 13. II In the same way, we obtain that 

\\T'V{^o + e-)\\HHB)<C\\V{^o + e-)\\mB), 

where C is independent of e. It remains to show that the operator / + e^T + e'^~'^^'^R 
is injective and to use the Fredholm alternative to obtain the existence of a unique 
w* G H^{B) verifying fl53l) . Injectivity is obvious when rj g]0, 2] since he operator norm 
of e'^T + e^+'^W in C{L^{B)) is of order 0{ei) < 1 for e < small enough. 

When ?7 = 0, we need to use assumption (H-2). Since —1 is not an eigenvalue 
of T, it suffices to fix Eq small enough such that the distance between —1 and the 
nearest eigenvalue of T is larger than £q~'^\\R\\c{l^{b))- To do so, we remark, following 
proposition 13.11 that uniformly in y G 5, -R(-,y) G W'^'^{B) with p < when 3 < 
d < 5 and p < oo when d = 2. The Sobolev embedding then yields that R{-, y) G C^{B) 
when d < 3 and R{-,y) G L'^{B) with q < oo when d = 4 and q = 5 when d = 5. Hence, 

||^1I/:(L2(B)) < C||i?(xo + £:-,Xo + £:-)IIl2, 

which is 0{1) for d < 3, 0{e~°') for any a > when d = 4, and 0{e~^) for d = 5. For 
the particular case go = 0, proposition 13. II gives R G C°°{B x B) so that ||-R^||£(l2(b)) is 
bounded independently of e for any d. In any event, eQ~'^\\R'^\\c(L2{B)) = o(eo) so that 
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the Fredholm alternative yields again a unique w* G L'^{B) solution to for small 
enough. In addition, w* satisfies the estimate: 

|k*||L2(i?) <C£''||y(xo + £-)llL2(B). (54) 

Then w'^ is given, for y G ^2, by: 



e 



(y) = [s^Tw* - e'^~^+'>R'w*) (^-7^) - ^'^(y). otherwise. 



so that w"^ G H^{Q). We verify that w'^ is then a solution to the variational formulation 
of ( l52l) . To prove uniqueness, we show that, for a given u G H^{Q), the assertion 

J S/u-\/vdx + J (^qo + ^q^(^-^^^ uvdx = 0, G if^(fi), (55) 

implies u = 0. Indeed, for ip G L^{^), consider the weak solution v G H^{Q) of 

—Av + qov = (p, X G n, 

augmented with homogeneous Neumann conditions on dQ. Thus, v is given by v{y) = 
Jj^ A^(x, y)v9(x)(ix. Plugging v into fl55|) leads to 



/" (m + T'u)ipdx = 0, G L2(1]), 



so that u + T'^u = 0, which implies that u = 0. This ends the proof of existence of a 
unique solution of the variational formulation of fl32]) when d > 3. 

We treat now the case d = 2. When 77 > 0, existence and uniqueness can be 
established in the same manner as above. When 77 = 0, we use assumption (H-3). We 
first notice that for d = 2, we have 

W(xo + £y)= / gi(x)^i;*(x)r(x-y)rfy-i^ / gi (x) u;*(x)rfx + i?W(y). 
Jb ^tt Jb 

In the same way, proposition 13.11 gives, uniformly in y G -B, i?(-,y) G W'^'''^'{B) C C^{B) 
with p < 00, so that we can recast R^w* as 



i?W(y) = i?(xo,xo + £y) / gi(x)ti;*(x)dx + i?W(y) 

Jb 

R'w*{y) = / gi (x) w*(x) (i?(xo + ex,xo + ey) - i?(xo,xo + £:y))(ix. 
Jb 

The system ( l53l) can then be reformulated as: 



+ Tu;* + R'w* = -TV - R'V - ( + i?(xo, xq + e-) ) C 
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where the constant is equal to 

= / Q'l (x) (^(xo + ex) + tf;*(x)) dx = / gi (x) i;^(xo + £:x)c/x. 
Jb J b 

Under assumption (H-3), we have = so that we just need to show that 

\W\\C{L^{B)) = 0{1) 

to apply the Fredholm alternative. Since i?(-,y) G C^{B), uniformly in y, we have, 
for all (x, y) E B x B, that |-R(xo + ex, xq + ey) — -R(xo,xo + ey)\ < Ce, which gives 
||-R^||£(L2{_B)) = ^'(e) and ends the proof of existence when d = 2. 

We now prove decomposition (124|) . which is the corner stone of the proof of theorem 
I3.3I Since w'^(xo + ex.) = w*{x.) when x G -B, it suffices to obtain an expression for w*. 
We consider first the case d > 3. Defining V^^(x) := ^(xo + ex), we recast (153|) as: 

w* + e'^Tw* = -e'^TV - e^-^+'^R' + w*). 

We expand in the first term of the right hand side and set w* = e^"^^ + e'^~'^^^ + ry, 
so as to obtain: 

e\J\ 



lil=o 

r' + e^Tr' = -R' {V + w*) , r'y + e'^Tr'y = -TR'y, 



where Ry is the remainder of the Taylor expansion of G C°°{B) of order d + 2. 
Writing vl/-(x) := Et^=ofd^Vi^o)<P]{^), with 

and following the preceding proof of existence when ci > 3, we verify that Vy G H^{B) 
with a norm bounded by Ce'^'^'^ and that and (jij are uniquely defined in H^{B). Also, 
examining ||-R^||£(l2(-5)) as in the proof of existence, we find that is bounded in L'^{B) 
independently of e when d = 3, is 0{e~°') for any a > when d = 4 and 0{s~^) when 
d = 5. When go = 0, is bounded in H^{B) independently of e since ||-R^||£(_ffi(B)) is 
uniformly bounded. We thus obtain the expression (!24l) announced in the proposition 
for d > 3. When d = 2, the equation for r*^ has to be replaced by 

r' + e'^Tr' = -R' {V + w*) + ^ / gi (x) t;^(xo + ex)dx, 

27r Jb 

= -R'{V' + w*)-(^^ + R{^o,^o + e-)^C', 

where is the same constant as before. When 77 > 0, we verify that G H^{B) with 
a norm of order loge. When t] = 0, assumption (H-3) implies = 0. Since i?^ is 0{e) 
in C{L'^{B), H^{B)), we deduce that is 0{e) in H^{B) since V is uniformly bounded 
in B and w* is bounded in L'^{B) according to □ 
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Proof of Theorem \3.'J[ We express in terms of V and the Green function A^, to 
obtain, a.e. in f2: 

^'(y) = ^(y) - e"-^^"" I gi(x) (V(xo + ex) + «;^(xo + ex)) iV(xo + ex, y)rfx. (56) 

Jb 

Taking the trace of (l56l) on which is well defined in LP'[dVi) and thus almost every- 
where, replacing tf^(xo + ex) by the expression in (12^ and Taylor expanding both V 
and N according to fl20|) . lead to the result. □ 



Proof of Proposition \3.4\ The outline of the proof is as follows: starting from the 
asymptotic expansion for v'^ in theorem 13.31 our aim is to recover that of in theorem 
12.21 and the expression of the polarization tensor M . This is done in several steps. First, 
we verify that assumptions (H-1), (H-2) and (H-3) are satisfied for the particular form 
(l26l) of the potential qi. In a second step, we show that the term in the expansion of 

is of order C(e^) so that 52(^-2) ye jg Q(^^2d'^ ^nd can be treated clS cL remainder. Then, 
we show in (l271)-(l28l) that the two first-order terms in the expansion of are actually 
of order 0{e^'^) so that they can be neglected and the expansions for and have 
the same leading order 0{e'^). Finally, using the particular form of the potential qi, we 
perform some transformations in the polarization tensors Q and for 77 = leading 
to the expression of the polarization tensor M in theorem 12. 2[ 

We will need the following lemma, which is one of the main ingredients to show the 
equivalence of the tensors: 

Lemma 4.1 Assume v G H^{B) verifies in the distribution sense, 



-Av + q,v = h mV'{B), gi(x) = ^^ o + ^^iW ^ h E L^B). (57) 

^/Do + Di(x) 

Then, for all (p G H^{B) and harmonic in B, we have: 

f gi(x)i;(x)<^(x)rfx=^ / Z}i(x)v(^=i^L=) ■ V<^(x)rfx - / /i^rfx. 



Proof Define L'(x) := L'o + ^i(x). Note that d^Di = on dB since Di G C^{n) 
and Di is supported in B. Hence, two successive integrations by parts yield: 



^^.(x)^(x)cix = - / vVD.v(^)d.= f (VD-^)A(^)d.. 



The above expression makes sense since (f is harmonic and Av G L'^{B) because of (l57l) . 
Starting from ( 1571) . we verify after some algebra that v solves 

V ■ Dvi^-^Jj =VDh, inV\B), (58) 



which, since D > in W^, is equivalent to: 



2V^D-V(^)+^a{^] = h. 
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Since D = Dq on dB and is constant, it follows from the above equation and another 
integration by parts that: 



■ V^^cix = / hifi(h(.. (59) 
J B 



Here, a is the surface measure on dB and the boundary term above has to be understood 
as the H~2{dB)-H'i{dB) duality product since d^v G H~2(^dB) because v G H^{B) 
and Av G L'^{B) thanks to fl571) . Using the fact that ip is harmonic in B, that D = Dq 
on dB, and using equation (jSHD, we find: 



B 



B 



D 
D 



C5 c/x + 2 



V(f (ix. 



B 



Vip (ix - 



V \ f dv 

■V(j9ax + 2 / -^—(pda 



h(pdx, 



dB on Jb 



'B WD 
To conclude, we just need to remark that, thanks to (158!) 



dv , 
— V^rfa 
dB dn 



1 



^0 JdB^dn 



(fda 



DV 



Jb 



V(fd:si. 



Coming back to the proof of proposition 13.41 we first verify that assumptions (H-1), 
(H-2) and (H-3) are satisfied. Since go = 0, (H-1) trivially holds because of the 
compatibility conditions ( ITTI) . The same is true for (H-3). Regarding (H-2), we have 
to show that if 

^ + Tip = 0, \/ipeL\B), (60) 

then ip = 0. To this aim, we first remark that T maps L'^{B) to H^{B), so that 
every (f verifying (!60|) belongs to H^{B). Now, ip can be extended to to a function 
if* G Hl^iW^) by the relation: 

^*(y)=v?(y), ye 5, 

(y9*(y) = —Tipiy), otherwise. 
Moreover, when (l60l) holds, then so does the following in the distributional sense: 

- Av?* + qi^* = 0, in P'(fi'), (61) 

for any bounded set Q' C W^. Consider y G M.'^\B. Then r(x — y) is harmonic for 
X G -B. We then apply lemma l4!T] with h = to find, uniformly in y: 



^*(y) = - / gi(x)v9*(x)r(x-y)dy, 

'B 



DA^W 



Jb 



Do + Di(x) 



Vr(x - y)rfx. 
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We thus deduce from the above equation for c/ > 2 the following behavior at infinity: 

¥^*(y) = 0(|yr-'), Vv.*(y) = 0(|yr^). (62) 
Besides, equation (I6T1) can be reformulated as: 



V- (Z^o + /^i)V 



Bji DD 5 of radius R and boundary Sr, we find: 

2 



0, in V\VL'). 



(63) 



After multiplication by (p*{Do + Di) 2 in iJj^^(]R'^), and an integration on the sphere 



{Do + Di] 



Br 



V 



5i, 



^(p*da = 0. 
an 



Letting R ^ 00 leads, together with fl62|) . to ^9 = so that assumption (H-2) is satisfied. 

We now show the equivalence of the tensors. First, the term given in the expansion 
of theorem 13.31 is of order O(e^), which is not obvious at first sight. Consequently, 
^2{d-2) j 1^^-^ is of order 0{e'^'^) and can treated as a remainder in the expansion. To prove 
this, we apply lemma 14.11 to and need to estimate r^. Let us recall the equation 
verified by G H^{B) given in proposition 13.21 



(y) + Tr^(y) = j qi (x) w^(xo + ex)i?(xo + ex, xq + ey)dy.. 

JB 



(64) 



When c? = 2, we use the fact that assumption (H-3) is satisfied since go = so that the 
term involving log e in the equation of proposition 13.21 vanishes. Since v'^ verifies (1571) 
with /i = 0, and R verifies (12T!) with go = so that we have Ayi?(x, y) = Ayi?(y, x) = 
since R is symmetric in its arguments and is thus harmonic, we apply lemma WA\ to find: 



/ gi (x)t;^(xo + £:x)_R(xo + £:x, xo + £:y)rfx 
Jb 

t;^(xo + £:x) 



e 

~ 7^0 J B 

Moreover, we show that 



/^i(x)V 



Vx-R(xo + ex, xo + ey)dyi. 



i;'^(xo + £■) 



0{e) 



(65) 



L^(B) 



SO that the left hand side of ( IMl) is of order 0{e^). This is obtained by proving that 
the leading term in the above expression vanishes. That is to say, thanks to the de- 
composition given theorem 13. 3[ w^(xo + ey) = V{'x.o + ey) + X]|^|io 7r^"'^(^o)0i(y) + 
£'i-V^(y) + C(£:'^+2), y a.e. in B, that 



^/ y(xo)(l + 0o(x)) x^ ^ 



(66) 



The argument is very similar to that in the verification of assumption (H-2) and so we 
just sketch the proof. Since 0o verifies 0o + = ~T1, it can be extended to to 
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\oc\ 



which admits the behavior at infinity given in (1^2]) . We also have, for any 



bounded set Vt' C W^, 

-A(0* + l)+gi(</.* + l) = O, inV'in') 
so that, still denoting by Br the sphere of radius R, 

(1)1 + I 



(67) 



Br 



V 



dyi 



Sr 



dn 



{ifl + l)da = 0. 



Sending R to infinity then gives the result thanks to the decay of V^g infinity. Owing 
to this result, the decomposition fl2^ . the fact that and belong to H^{B), and 
is at least an 0{e) when ci = 2 as mentioned in theorem (13.31) . we get that (jHSD holds. 
Furthermore, using again the fact that R is harmonic, we verify from that solves 
in the distribution sense: 

-Ar' + qir" = 0, mV'{B). 

We cannot apply lemma W7U directly to flM|) since for (x, y) ^ B x B, we have 

-A,r(x-y) = 5(x-y), mV'{B), 

and r is not harmonic. Nevertheless, the lemma can easily be adapted to this special 
case so that, y a.e. in B, we have 



gi(x) r^(x)r(x - y)dx 



1 



B 



"Do 



B 



Vr(x-y)dx 



VDo + A(y)-v^ 



Plugging the above expression into (IMI) . we finally find the following equation for G 
H^iB), y a.e. in B: 



+ — / Z}i(x)V 



/^fl.(x)v( 



Vr(x-y)c/x 



^Do + /^i(x) 

• Vx-R(xo + ex, xq + ey)dx. 



+ A(x) 

Identifying the right hand side of the latter equation with S^{ex) and (i5o+-Di)~5r^ with 
(x) + S'^(£:x) in the proof of theorem [221 we see that (Do + -Di)~2r^ and rf (x) + 5'^(£:x) 
satisfy similar equations so that the same technique yield 



<Ce' 



LHB) 



LHB) 



I VxVyi?||^oo(5QXBo) 



From proposition 13.11 R G C°"{Q x fi). Together with ( l65|) . this finally gives that: 



V 



(68) 



L2(B) 
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We conclude by applying once again lemma 14.11 to obtain 



/ gi(x)r^(x)iV(xo + ex, ■)dx. 
Jb 



/Dn 



/^i(x)V, ^ 



VxiV(xo + £x, ■)rfx 



L2{sn) 



thanks to (l68l) . 

We now prove (127!) and (128!) so that the leading order in the expansion of theorem 
12.21 is 0{e'^) as in the case of the diffusion equation. We remark that, for rj = 0, 



= gi(x) (^^(x) + y(xo + ex) - i?^v(x))dx, 

where is given in the theorem and Ry is the remainder of the Taylor expansion of 
V^(xo + ex) at the order d + 2 and is thus of order 0{e'^^'^). In order to apply lemma 
I4n we verify from ([25]) that G H^{B) solves, 

-A*^ + gi^^ = -gi (v(xo + ex) - /2^(x)) in V'{B). 



Setting v{-x.) = \l'^(x) + V"(xo + ex), h = qiRy and (^9 = 1 in lemma 14.11 yields (1271) . 
Regarding (!28|) . we write, for y G dfl, 

d+l |j| „ 

^i_a»7V(xo,y)(Q,o + Q°o) = y^gi(x)(l + 0o(x))(iV(xo + ex,y)-i?^(x,y))rfx, 

where R^j is the remainder of the d + 2 order Taylor expansion of A^(xo + ex, y) with 
respect to x and is thus of order 0{e'^~^'^). Since A^(xo + ex, y) is harmonic when x G -B 
and y G dQ, we apply lemma 14.11 thanks to (!57|) to find: 



J]^aw(xo,y)(g.o + g°o) 

1=0 ^' 



;i+0o(x)) 



+C)(e''+2) = C(e^+2^, 



VxiV(xo + ex,y)c?x 



since the above integral vanishes thanks to 

At this point of the proof, we have thus shown that satisfies, a.e. on dQ, that 



d+l d+l £;_2+|j|+|j| 



2d\ 



|j|=l|i|=l 



Setting f^(y) := n^(y)y'-Do + -Di(^^-^), V := y/D^U, we verify that and U are 
solutions to ([T]) and (jlj), respectively, with the boundary term g multiplied by y/D^. It 
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thus remains to show that fl2U]) and fl30p hold to recover the asymptotic expansion for 
u"^ of theorem 12. 2[ Since x-' + (pj satisfies fl57j) when | j | = 1 and x* is harmonic when 
\i\ = 1, we have, for \i\ = \j\ = 1: 



Qij + Qi'. 



gi(x) (x.^ + (pji^j^iddx, 

xJ + 0,-(x; 







^Z}o + ^i(x) 



Vx'rfx. 



(69) 



We introduce the following extension to 0,- on 



ye 5, 



fjiy) = -T(t)jiy) - Tx.\ otherwise, 

which thus satisfies the conditions at infinity in (1621) . We recall that (p^Q, the function 
introduced in theorem 12.21 to define the polarization tensor M, is the unique weak 
solution in the space Hl^^(R'^) fl C°°(M'^\S) to the following system posed in R*^: 



-V- (^Di(x)Vx^ 
0(|x|^-^ 



as X 



oo. 



(70) 
(71) 



When \j\ = 1, notice that (f)% is given by 



- 1 X' 



x), 



so that ( j29l) is proved using (j69l) . To prove (j30l) . we need to sum over i and j to be able 
to use lemma HIT] since x* is not harmonic for |«| > 2 and x-' + 0j(x) satisfies ( |571) with 
a negligible left- hand side h of order (^(e'^"''^) only after summation. We thus write, 
using the same arguments as for the proof of fl27|) and (l28l) . for y G 

\j\=i\i\=i 

= EE 

iii=iiii=i 



9W(xo,y)a^V(xo) / gi(x) x^ + 0,(x) x^rfx 



/ gi(x) (y(xo + ex) - i?^(x) + ^^(x)) (iV(xo + ex, y) - i?^(x, y))(ix + 0{e''+^), 

JB 

■■l/(xo + £X) + ^^(x) 



Z^iV 







Z}o + ^i(x) 



■ VxiV(xo + ex, y)rfx + 0{e 



d+l d+l + 



EE 

b1=i Nl=i 



-9W(xo,y)9^V(xo) / A(x)V 



x^ + 



■Vx^dx+ C(e 



It remains to relate the latter sum to M. For that, let fj be defined as: 

My) 



<P*iy)-<i^l{y)- 
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Then fj belongs to Hl-^^lW^) fl C°°(M.'^\B) by construction and solves: 

V-{Do + D,{^))Vf, = -]Ib(x)v//^o + ^i(x) Ax^ x G (72) 
/,(y) = Oi\y\'~') as |y|-oo. (73) 

Here, Us is the characteristic function of the set B and 0* is the extension of (pj to M'^. 
Note that fj = when | j | = 1 so that we recover the preceding relationship between 0* 
and 0Q. To conclude the proof, it suffices to show that an appropriate linear combination 
of the terms fj is of order (9(e'^+^). Let: 

d+l |j| d+1 |j| 

^v(x) := E '-rd^Vi^o)A^', i^^(x) := ^^^■^(xo)/.(x), 
lil=i lil=i 

so that since AV{'Kq + ex) = 0, for all x G 5, we have Ty(x) = 0{e'^~^'^) uniformly in B 
and G i/ioc(K'^) H C°°(M'^\E) solves 

V-(Z^o + ^i(x))VF^ = -]Ib(x)v/^o + A(x)T^, xgR^ 
F%y) = Oilyf") as |y| ^ oo. 

The above equation is very similar to fj45l) at the end of proof of proposition 12.81 and a 
similar analysis yields 

||VF^|U.(M.) = 0(e'^+'). 
We conclude the proof by calculating that 

d+1 d+1 \i\ + \j\ f 3 A_ A ( \ 

d+1 d+1 \i\ + \j\ 1 f 

= E E ^^^A^(xo,y)^V^(xo)^ / /^i(x)v(x^- + + /,) ■ Vx^cix, 
1,1=1 |.|=i ^^^Jb V / 

= ^ E E -rr^^^(^0'y)^'^(^o)^^^- + )■ 

□ 

4.3 Appendix 

This appendix states several lemmas that were needed in the preceding analyses. 

Lemma 4.2 Let F G {L'^{W^)Y and Di G W^'°^{W^) compactly supported in a hounded 
domain B, and Do a strictly positive constant. Assume moreover that Dq + -Di(x) > 
Co > a.e. in M.'^. Then, the following problem (P): 

V- (1^0 + A(x)) V0 = V-F mV'iW^), 

0(x) = 0{\^\^-'^) as |x| ^ oo. 
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admits unique solution in H]^J^'^)r\C'^{W\B). Moreover, (f) satisfies the estimates, for 
any hounded set A C , 

||V0||L2(igd) < Co"^||F||(i2(B))d, ||<^||l2(a) < C'l|F||(L2(B))d (l + ||Di||ioo(Kd)) , (74) 

and is the unique solution, a.e. on every hounded set ofW^, to the integral equation 



= - Di(x)V</.(x)- Vxr(x-y)rfx+ / F(x)- Vxr(x-y)dx. (75) 

J B J B 

Proof. We show that (P) is equivalent to a problem posed on a bounded domain 
that can be solved with the Lax-Milgram lemma. To do so, let be the sphere of 
radius R with B GG Br and denote by Sr its boundary. Consider the solution to (P) 
with the announced regularity. Since both Di and F are supported in the function 
is harmonic in W^\B and in particular in W^\Br. Denoting by A : H2{Sb) H~'2{Sb) 
the exterior Dirichlet-Neumann operator on the sphere Sr, we then have the standard 
relation 

where |^ is the outer normal derivative of on Sr and ^l^.^ its outer trace. Since is 
harmonic in W''\B and is thus of class on this set, |^ and 0|^^ are continuous across 
Sr. Using this fact and integrating (P) against a test function v G C°°{Br), we find 

/ {Do + D,)V(j)-Vvd^-Do{A(j)\s^,v\sJ= [ F-Vwdx, 
Jbr Jb 

where (■, ■) denotes the H^{Sr) — H~^{Sr) duality product. The restriction of to Br 
is therefore a solution to the following variational problem (P2): Find u G H^{Br) such 
that 

a{u,v)=l{v), yveH\BR), 

with obvious notation for the bilinear form a and the linear form /. Let us assume for 
the moment the existence of a unique solution u to (P2). That solution can be extended 
to a function u* solution to (P) . Let indeed u* be defined as: 



u* = u, in Br, 



u* = U, in R'^\Bt 



where U is the solution to 



AU = in P'(R'^\Ph), 
U\sr ~ '^^Sr^ t/(x) — > as |x| — * oo. 

By construction, the trace of u* is continuous across Sr. Since U is harmonic in M.'^\Br 
and vanishes at infinity, it also verifies: |^ = AU\g^ = Au\g^ . It then suffices to 
integrate the equation solved by U against a test function v G C^(M'') and to consider 
(P2) to find 



/ {Do + Di)Vu* ■Vvd^= I ¥-Vv c/x, G C, 

jRd J B 
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so that u* solves (P). The above equation also implies that u* is harmonic in W^\B 
and is thus of class C°° on this set. It remains to verify the behavior at the infinity, 
which stems from the fact that F has compact support in Sr. Setting w = 1 in (P2) 
yields (A , 1) = 0. Getting back to U, since its trace and its normal derivative are 
known and given by u\g^ and Au\g^, respectively, it admits the following representation 
formula, for x G M.'^\Bji: 

f/(x)= / u\s, (y) ^^^^ ~ ^^ da{y) - (A ^|^^ , r(x - ■)), 

JSr oily 

where F is the fundamental solution of the Laplacian in and a is the surface measure 
on Sr. We conclude by noticing that, as |x| — > oo: 



(A , r(x - •)) = (A u\,^ , r(x - ■) - r(x)) = 0(|x 



l-d\ 



It remains to show the existence of a unique solution to (P2). This is a consequence 
of the Lax-Milgram lemma: a and / are both continuous in H^{Br) and the coercivity 
follows from the Poincare-type inequality: 

Mlhbr) < C {\\Vu\\l2^Br) + Mlhsr)) , Vn G H\Br), 
and the relation 

We now prove the first estimate in (CI). Let G C^{R'^) such that 

■^MIlhsu) (76) 
Integrating (P) against v yields 

[ {Do + Di)\/(f)-Wvd^- Do [ ^vda= [ F-Wvdx. 
J Br Jsr dn Jb 

Since V</)(x) = Cdxl^"^) as x tends to infinity, it belongs to L'P{W^\Bp) for some p > \ 
and a ball of radius p with B GG Bp. The above equality also holds by density for all v G 
Vp, the space of functions v such that v G Hl^^{R'^), v verifies ( !76|) and Vv G (M.'^\Bp) 
for ^ + ^ = 1. Since |^ = 0{R~'^), sending R to infinity implies, together with (j76ll . that 
the boundary term goes to zero. On the other hand, the function V0 ■ Vw is integrable 
on M*^ for v G Vp, which allows us to use the Lebesgue dominated convergence theorem 
and obtain as i? — > oo: 

[ {Do + Di)V<j)-Vvd^= I F-Vvd^, (77) 

JlRd Jb 



for all V G Vp. Since G for any o? > 2, we obtain the left estimate of (174|) . 

Let us now consider the integral equation (!75|l and show that the solution to (P) 
verifies ([75]). For G L'^{Br), let w(x) = F(x-y)^/;(y)dy for a given ball Br. Since 

F G Wl^l{W^), it follows from the Young inequahty that v G Hl^{W^). Set x G W^\B^ 
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with Br CC Br,. Then Vr(- - y) G Lp{W^\Br,) for p > ^ and y G Br, Such 
a function v also satisfies fl76l) for d > 2 since r(x — y) grows at worst as log |x| for 
(x, y) e W^\Bri X Br. We can thus use w as a test function in fl77|l . In order to use the 
Fubini theorem, we notice that the function Vt'(x) ■ Vr(x — y)-?/'(y)IlB^(y) belongs to 
L^(M'^ X W^) thanks to the Sobolev inequality [TT] recalled in lemma in the appendix 
since Vw G L^(M"') and ip G LP'^Br). Indeed, since R < oo, we bound the L^(M"') norm 
of '?/^(y)IlB^(y) by the L'^{Br) norm of ip for q = < 2. Then choose p = 2 and 
A = (i — 1 in lemma 14.31 

The same conclusion holds for F(x) ■ Vr(x — y)'?/;(y)IIs^(y) so that we obtain from 

Do [ (I V0(x). vr(x-y)dx)v(y)c^y 

J Br ^Jk'' ^ 

Di(x)V0(x) ■ Vr(x - y) - F(x) ■ Vr(x - y))iP{y)dy. (78) 



'B 



R 



It thus only remains to show that J^^Vcplx.) ■ Vr(x — y) rfx = </'(y) a.e. on Br to 
conclude. To this aim, consider a sequence (/)"■ of [W^) functions such that V^" — > V0 
in L^(M'^) and 0" ^ in L'^{A) for any bounded set A. Since — Axr(x — y) = 5(x — y) 
in the distribution sense, we have, for any y G M'^: 

lim / r(x - y)A0"(x)dx = -0"(y). 

The Lebesgue dominated convergence theorem yields consequently: 

hm [ ( [ r(x - y)A0"(x)rfx') ^(y)rfy = - / 0"(y)z^(y)rfy. 

An integration by parts then gives: 



[ r(x-y)A0"(x)dx = [ ^^r(x-y)rfa(x) 

- [ Vr(x-y) ■ V0'^(x)c/x. 



'|x-y|>£ 

The boundary integral goes to zero with e. For the other term, we remark that the 
function ]I|x_y|>elIijBVr(x - y) ■ V</)"(x)?/'(y) converges a.e. in R'^ x M'^ to ]Iij^Vr(x- 
y)- V</'"(x)^/^(y) which belongs to L-'^(]R'^xM'^) thanks to the Sobolev inequality. Applying 
again the Lebesgue dominated convergence theorem yields 

[ ( [ V0"(x).Vr(x-y)rfx')^(y)rfy= / 0"(y)V^(y)c^y, 

JBr XJRd. / J Br 

and it suffices to pass to the limit in the sequence 0" to conclude. This proves that the 
solution to (P) satisfies fl75|) . Conversely, considering a solution of fl75|) in Hl^^iW'-) fl 
C°°(M'^\i?), we verify using the same techniques as above that this solution also satisfies 
(P), which we know admits a unique solution. Therefore, the integral equation fl75l) also 
admits a unique solution. The second estimate of fl74|) follows from flTSl) . the Young 
inequality and the first estimate of fl74p . □ 
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Lemma 4.3 Sobolev inequality (see e.g. 177]/). Let f G Lp(R"'), g G L''(R'^), 1 <p,q < 
oo, < X < d with the relation - + - + 4 = 2. Then: 



The following lemma, which is a standard variational formulation of the Fredholm al- 
ternative, is used several times in the paper. 

Lemma 4.4 Let H be a Hilbert space and let a{-, ■) be a bilinear form on a H x H 
such that a(-, ■) = ao(-, ■) + ai(-, ■), where both and ai are continuous in H and is 
H -coercive. Assume moreover, that for two sequences Un and w„ weakly converging in 
H to u and v, we have 

ai{un,Vn) ai{u,v). 
Then, if the following assertion is verified 

{a{u, v) = 0, e H) ^ u = 0, 

for all f in H' , there exists a unique u E H which satisfies 

aiu,v) = {f,v), \/veH. 

Here, (■, ■) denotes the H'-H duality product. Moreover, u verifies the estimate, for 
some positive constant C: 

\\u\\H<C\\f\\H'. 

Proof. We sketch a proof for completeness. Since Oq is coercive, we know from the 
Lax-Milgram theory the existence of a bounded and boundedly invertible operator S 
on H such that ao(u, t') = {S~^u,v), where (■, ■) is the inner product on H. By the 
Riesz representation theorem, we similarly know the existence of a bounded operator 
Ai such that ai{u^v) = {Aiu,v). The hypotheses on oi imply that Ai is compact on 
H. Indeed, choose Un ^ u and define Vn = AiUn — Aiu. We verify that f „ ^ and 
that II Ai^n — AiMp = {AlUn, Vn) — {Aiu, Vn) couverges to by the above hypothesis on 
ai so that Ai maps weakly converging sequences to strongly converging sequences and 
is thus compact. 

Now by the Riesz representation theorem, there exists f E H such that (/, v) = 
{f,v), for all V E H, so that a{u,v) = {f,v) is equivalent to [S^^ + Ai)u = f and 
thus equivalent to (/ + SAi)u = Sf, which admits a unique solution if and only if —1 
is not an eigenvalue of the compact operator SAi, which is equivalent to the fact that 
a{u, v) = for all v E H implies that u = 0. □ 
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